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ABSTRACT 

We have determined the luminosity function of 250yum-selected galaxies detected in the ~14deg^ science demonstration region of the Herschel- 
ATLAS project out to a redshift of z = 0.5. Our findings very clearly show that the luminosity function evolves steadily out to this redshift. By 
selecting a sub-group of sources within a fixed luminosity interval where incompleteness effects are minimal, we have measured a smooth increase 
in the comoving 250^m luminosity density out to z = 0.2 where it is 3.6^Qg times higher than the local value. 
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1. Introduction 

■ Measurement of the galaxy luminosity function (LF) constitutes 
, one of the most fundamental statistical constraints that can be 

■ placed on models of galaxy formation, and hence the build up 
of large scale structure, in the universe. Since half of the en- 

, ergy ever emitted by galaxies has been absorbed by dust and re- 
' radiated at far-infrared and sub-millimetre (submm) wavelengths 
' (Fixsen et al. 1 19981 1 and, because knowledge of the statistical 
properties of submm sources is relatively sparse, determination 
of the submm LF provides a crucial missing piece in a fully com- 
I prehensive model of galaxy evolution. 

Following their detection in the first deep submm and mm 
\ surveys (e.g.. Small et al. 1997; Hughes et al. 119981 Eales et 
. al. [1999; Bertoldi et al. 12000k much has been learned about 
' the dusty high-redshift sources selected at such wavelengths. 
\ Although many studies have argued that these sources are likely 
I ancestors of local ellipticals (e.g., Scott et al. 120021 Dunne et al. 
[2003 ) little progress in verifying this assertion has been made 
since. The main reason for this is the preponderance of high red- 
shift submm/mm-selected sources, owing to the strong negative 
k-correction and small survey areas. The resulting low numbers 
of sources at redshifts z < 1 has therefore precluded evolutionary 
studies over the last ~60% of the Universe's history. 

In particular, despite the local submm LF being first deter- 
mined a decade ago (Dunne et al. 12000 J . little has been added 
to our comprehension of how the LF has evolved over the 



last ~ 7 Gyr, until very recently. Observations conducted using 
the Balloon-borne Large Aperture Submm Telescope (BLAST; 
Devlin et al. '2009") have made significant improvements with 
much enhanced sensitivity to z < 1 sources. As a result, direct 
estimates of the LF at 250, 350 and 500 //m were made by Eales 
et al. (120091) who detected strong evolution, particularly among 
the higher luminosity systems, from z = 1 to the present day. 
However, the accuracy of these findings is limited by the small 
number of sources used (~ 50 at z < 0.5) and source confusion 
due to the angular resolution of BLAST. 

In this letter, we present our measurement of the LF 
of 250yL(m-selected galaxies detected by the Herschel Space 
Observatory (Herschel; Pilbratt et al. 2010) over a ~ 14deg^ re- 
gion acquired as part of the science demonstration observations 
of the HerscheJ-Astrophysical Terahertz Large Area Survey (H- 
ATLAS; Eales et al. 120101) . These data offer a significant im- 
provement over the BLAST data in terms of their increased sen- 
sitivity, higher angular resolution and greater areal coverage, re- 
sulting in ~ 20 times the number of sources with which to com- 
pute the LF. 

Throughout this letter, the following cosmological parame- 



ters have been assumed; 
Qa = 0.73. 



2. Data 



Ho 



71kms"'Mpc-', Qn 



0.27, 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA 



The 4° X 4° H- ATLAS science demonstration field was observed 
with the Spectral and Photometric Imaging Receiver (SPIRE; 
Griffin et al. 2010) at the wavelengths 250, 350 and 500 //m and 
with the Photodetector Array Camera and Spectrometer (PACS; 
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Poglitsch et al. l20T0b at 100 and 160/im. The field, centred at 
the co-ordinates (09''05"'30', +00°30'00"), was scanned twice in 
parallel mode. The 5cr point source sensitivities of the resulting 
beam-convolved maps, including confusion noise, are 132, 126, 
32, 36 and 45 mJy and the beam sizes expressed as full width at 
half maximum (FWHM) are 9", 13", 18", 25" and 35" at 100, 
160, 250, 350 and 500 yum respectively. Details of the SPIRE and 
PACS map-making are given in Pascale et al. (in preparation) 
and Ibar et al. (in preparation) respectively. 

Sources were initially extracted from the 250 /^m noise- 
weighted beam-convolved map in a central 14.4 deg^ region 
above a significance of 2.5cr. For each source, 350 and 500 //m 
fluxes were then estimated from the appropriate beam-convolved 
map at positions determined at 250fj.m. Extended source fluxes 
were measured in apertures matched to identified (IDed) optical 
counterpart sizes (see below). 6878 sources were detected with a 
significance of > 5cr in any one band. Fluxes at 100 and 160yum 
were assigned by matching to > 3cr PACS sources within a posi- 
tional tolerance of 10". The analysis described hereafter applies 
to the 6613 sources detected at > 5cr at 250 ;um. Full details of 
the source extraction are given in Rigby et al. (in preparation). 

Since our field (and H-ATLAS at large) is lacking in the ra- 
dio and mid-infrared data traditionally used to identify counter- 
parts to submm sources, we have taken a different approach and 
matched directly to optical counterparts. It is possible to attain a 
reasonable rate of secure optical IDs in this way with H-ATLAS 
because, unlike existing submm surveys at 850yL(m-l.l mm sen- 
sitive to high redshift sources (Zmedian - 2.5), H-ATLAS sources 
lie at substantially lower redshifts (Zmedian $ 1 -0) on average. This 
means that we can use shallow optical imaging, where a low sur- 
face number density of sources allows for much less ambiguous 
IDs. 

We used the likelihood ratio (LR) method of Sutherland & 
Saunders ( II992I I to perform the matching, which uses the submm 
positional uncertainties and the magnitude distribution of coun- 
terparts to assign the likelihood that a particular optical source is 
physically associated with a target submm galaxy. We searched 
for optical counterparts at r < 22.4 from the Sloan Digital 
Sky Survey (SDSS) seventh data release (Abazajian et al. 2009) 
within 10" of every 250 ;um SPIRE source. The LR technique as- 
signs a reliability parameter, Rlr, to each match, which indicates 
the probability that the counterpart is the correct ID. The calcu- 
lation of Rlr includes the probability that the true counterpart 
may be below the detection limit of the survey and accounts for 
other counterparts within the same search radius. To remove un- 
reliable counterparts, all those with Rlr < 80% were discarded, 
leaving a total of 2267 submm sources with unique optical coun- 
terparts. We refer the reader to Smith et al. (in preparation) for 
an exhaustive account of the ID procedure. 

Of these 2267 counterparts, 876 have spectroscopic red- 
shifts (spec-zs) acquired either by the SDSS, the Galaxy And 
Mass Assembly survey (Driver et al. I2009I I. the 2dF redshift 
survey (Colless et al. |2001) or the 6dF redshift survey (Jones 
et al. I2009I I. For the remaining counterparts, photometric red- 
shifts (photo-zs) were estimated by applying the ANNz neural 
network code (Collister & Lahav '2004) to the SDSS optical 
photometry and also near-infrared photometry taken from the 
seventh data release of the UKIRT Infrared Deep Sky Survey 
(UKIDSS; Lawrence et al. 2007). Details of these photo-zs are 
given in Smith et al. (in preparation). A total of 2239 sources 
were assigned photo-zs, with a total of 2241 sources having a 
redshift of either type. Figure[T]shows the redshift distribution of 
these 2241 sources. The 1688 sources at z < 0.5 form the sample 
to which we apply our analysis in this letter In every case, we 
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Fig. 1. Redshift distribution of the 2241 sources that are detected 
at 250 /vm with > 5cr significance, and have an optical counter- 
part with a redshift and a reliability of association of > 80%. The 
1688 of these at z < 0.5 form the sample analysed in this letter. 

used a spec-z in preference to a photo-z, although there is ex- 
cellent agreement between the two, with a standard deviation of 
(Zphot - Zspec)/(1 + Zspec) = 0.039 ovcr the sample. 

3. The luminosity function 

3.1. SED fitting 

A modified black-body spectral energy distribution (SED) was 
fitted to the SPIRE photometry, and, where available (272 
sources), PACS photometry, for each source. In the fitting, we 
allowed the rest-frame dust temperature, T, to vary between 
lOK < T < 50K, we fixed the dust emissivity index to /? = 1.5 
and we fixed the redshift to either the photo-z, or, preferentially, 
when available, the spec-z. The temperature of 342 sources 
(none with PACS photometry) could not be reliably constrained. 
For these, we re-fitted the SED, fixing the temperature to the 
median of the sample (see below). 

Excluding the sources where T was fixed, we found a me- 
dian dust temperature across the whole sample of 26K with a 
standard deviation of 8K (or 23 + 7 K for fi - 2). This result 
closely agrees with the rest-frame dust temperatures determined 
for BLAST sources measured by Dye et al. (2009) as well as 
those of the sample of H-ATLAS galaxies studied in Amblard et 
al. (12010: see also other temperature comparisons therein). 

For each source, we measured the rest-frame 250 jum lumi- 
nosity, L250, by integrating the rest- frame SED over the 250 jum 
SPIRE bandpass function. Errors on L250 were determined by 
propagating redshift errors and temperature errors obtained from 
the SED fit. 



3.2. Estimating fhe LF 

We based our measurement of the LF on the estimator 



<f> = 2:,(ymax,i) 



(1) 



where Vmax.i is the comoving volume out to the maximum red- 
shift that source / could be placed and remain above the 250 /vm 
anc/ optical r-band detection threshold. We used an Sb type SED 
for the optical k-correction and the SED fit as described in the 
previous section for the k-correction at 250 //m. The sum here 
acts over all sources in a given luminosity and redshift bin. We 
computed the total error on as the quadrature sum of its for- 
mal error, ^J%iVmmJ^^, and the standard deviation of the scatter 
measured in each bin in performing a Monte Carlo simulation in 
which redshifts and photometry were randomised according to 
their errors. 
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then calculated as the reciprocal of the fraction of counterparts 
that would be IDed within an offset of the smaller of r^o and the 
ID search radius of 10", i.e.. 



100 

flux / mjy 



1000 



Fig. 2. Rate of identification of optical counterparts to 250 ;um 
sources (dashed line) and ratio of the number of z > 0.5 sources 
to the number of z < 0.2 sources (continuous line), as a func- 
tion of 250 jum flux. Both quantities were determined using the 
sample of 2241 counterparts. Poisson errors are plotted. 



An important consideration when computing the LF is in- 
completeness. There are two types of incompleteness at play in 
our case. Firstly, 250 ;um sources with optical fluxes lower than 
the r-band detection threshold are missing. This is a stronger 
effect at higher redshifts and lower luminosities. Figure|2]shows 
how the ID rate (i.e., identifying with an optical counterpart with 
Rlr > 80%) varies with 250 yum flux over the full sample of 
6613 sources detected at 250 ;um. The figure also shows that the 
ratio of the number of z > 0.5 sources to the number of z < 0.2 
source^ increases rapidly towards low 250jum fluxes. Together, 
these two facts imply that the unidentified 250 //m sources pre- 
dominantly lie at high redshifts where their optical flux falls be- 
neath the sensitivity limit of the counterpart catalogue. 

However, the fraction of missing sources at z < 0.5 (where 
we have computed the LF - see Section l4Tb with undetected op- 
tical flux must still be quantified and accounted for. We therefore 
plotted the optical r-band flux against 250 yum flux for sources in 
different redshift slices and found a clear positive correlation in 
every slice. At redshifts z > 0.2, the faint end of the locus of 
points on this plot becomes cut by the optical sensitivity limit. 
Using this result, we were able to determine a 250 jum flux limit 
for each redshift slice where we estimate that < 5% of sources 
are missing due to the optical sensitivity. We then limited our 
computation of the LF in each redshift bin by the correspond- 
ing 250 jum luminosity limit, ensuring > 95% completeness at 
all redshifts. H-ATLAS completeness effects will be discussed 
at length in forthcoming work, including Rigby et al. (in prepa- 
ration). 

The second incompleteness effect is a consequence of the re- 
liability cut. Approximately one half of the 250 yum sources that 
were matched to an optical source, were rejected by applying a 
minimum reliability of 80%. A certain fraction of these will be 
genuine counterparts. To estimate this fraction, we applied a cor- 
rection factor This factor was computed by reversing the algo- 
rithm used in the ID procedure (see Smith et al., in preparation) 
that determines a counterpart's reliability from the radial off- 
set between the submm and optical positions, the optical r-band 
magnitude and the signal-to-noise ratio of the submm source. We 
therefore computed the submm-optical offset, rgo, correspond- 
ing to a reliability of 80% for the full range of combinations of 
250 yum source signal-to-noise ratio and r-band magnitude seen 
in the data. For each combination, the correction factor, c, was 



-min(rgo,10") / ^2 \ 

c - A[ I rexp ^ dr 



Jo 






(2) 



' These two redshift limits were chosen to obtain approximately 
equally sized sub-samples. 



where the normalisation A is set such that one would obtain c - I 
were the integral evaluated between and oo. Here, cr depends 
on the submm signal-to-noise, fi, and the beam FWHM accord- 
ing to cr = 0.6 * FWHM/yU, as given by Ivison et al. ilOOIi . At 
250 /im, cr therefore varies from ~ 2.2" for a source with /i = 5 
to ~ 1.1" for a source with yu = 10. These values are consis- 
tent with the distribution of offsets obtained in matching to the 
SDSS. The minimum value of cr was limited to 1 .0", to account 
for the SPIRE pointing error and map pixel size (see Smith et 
al., in preparation, for more details). 

The resulting correction factor was then applied by modify- 
ing the LF estimator to 2,- c,(ytnax,i) ' where c,- is the correction 
factor corresponding to the ith source's 250 /im signal-to-noise 
and r band counterpart magnitude. The average correction factor, 
weighted by the counterpart number counts, ranges from ~ 1 .2 
for 5cr 250 /^m sources, through ~ 1.1 for 8cr sources to ~ 1.0 
for 1 Oct sources. 

Finally, we estimated the expected number of false counter- 
parts by summing the quantity 1 -Rlr- Within the 1688 sources 
at z < 0.5 used for computation of our LF, we estimate a total 
of 81 false counterparts. This false ID rate of ~ 5% shows no 
noticable correlation with redshift or luminosity and little vari- 
ation between LF bins across the redshift-luminosity plane. We 
treated the false counterpart rate as an additional source of error 
and added it in quadrature to the error on each LF bin. 

4. Results 

4.1. LF evolution 

Using the methods outlined previously, we determined the LF in 
five redshift bins of width Az = 0.1 from z = to z = 0.5. The 
total error budget includes the formal error on the LF estimator 
that accounts for Poisson noise, false IDs and the scatter mea- 
sured in the Monte Carlo simulation which randomises redshift 
and photometry. 

Figure [3] shows the LF in each of the five redshift bins. The 
figure clearly shows that the LF exhibits significant evolution out 
to z = 0.5. At a given luminosity, the comoving space density in- 
creases steadily with redshift. This is consistent with the findings 
of Eales et al. (2009), although our detection of evolution is con- 
siderably more significant. The figure also shows that the local 
250jum LF predicted by Serjeant & Harrison ( I2005I I agrees very 
well with our z < 0. 1 LF. 

In each redshift bin, the slope of the LF becomes shallower 
toward lower luminosities. Our allowed incompleteness of up to 
5% is insufficient to account for the magnitude of this effect, 
although it is possible that there may be additional mild incom- 
pleteness in the source extraction process at low 250 yum fluxes. 
This will be quantified in a later paper once the source extraction 
has been formally characterised (Rigby et al., in preparation). 



4.2. Luminosity density evolution 

Having established significant evolution of the LF, we investi- 
gated evolution of the comoving luminosity density. Clearly, in- 
completeness will preclude an accurate measurement out to any 
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Fig. 3. The 250 fim LF in five different redshift bins of width Az = 0.1 from z = to z = 0.5. Errors account for Poisson noise, 
redshift and photometric errors and the expected false ID rate. The thicker grey line shows the local 250 //m LF predicted by Serjeant 
& Harrison Q2005I I. 
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Fig. 4. Evolution of comoving rest-frame 250 ;um luminosity 
density out to z ^ 0.2 for sources with luminosity in the range 
10^ < L250/L0 < 5 X 10^. Errors were determined using a Monte 
Carlo simulation which randomised redshifts and photometry 
and account for Poisson noise. The line is the fit (1 + z)^ '. 



5. Conclusion 

One of the key goals of H-ATLAS will be to understand the na- 
ture of the evolution detected in this letter In turn, we aim to im- 
prove our understanding of the evolutionary link between high 
redshift and local submm systems. 

This letter has only considered sources selected at 250 yum, 
merely one of the five wavebands on offer from H-ATLAS. 
Furthermore, the 14deg^ of survey data analysed in this work 
represent only 2.5% of the final, proposed H-ATLAS survey 
area. A repeat of the analysis presented here with the final survey 
data, would therefore result in the quoted uncertainties falling 
by at least a factor of five. In light of these considerations, it is 
clear that our results offer only a small glimpse of the anticipated 
wealth of science that H-ATLAS has to offer. 

Acknowledgements. SD Acknowledges the UK STFC for support. 



significant redshift. However, by limiting the calculation to a 
sub-group of sources within a fixed luminosity interval where 
incompleteness is small, it is possible to estimate the strength of 
evolution in the sub-group over a larger redshift range. 

We therefore computed the comoving rest-frame 250 yum lu- 
minosity density of a sub-group of sources with luminosities in 
the range 10^ < L250/L0 < 5 x 10^, allowing measurement up to 
z - 0.2 before incompleteness becomes significant. We note that 
within this luminosity and redshift range, approximately 85% 
of the sources have spectroscopic redshifts with negligible er- 
rors. Figure ID shows the results. The luminosity density exhibits 
a steady and significant increase with redshift, which when fit 
with the form (1 + z)" yields a value of n = T.l^j^. This corre- 
sponds to an increase in luminosity density by a factor of 3.6^,, g 
from the present day to a redshift of z = 0.2. This scaling is 
consistent with Saunders et al. (11990 ^) who measured a density 



scaling of (1 + z) 



6.7±2.3 



out to z = 0.25 for 60jum-selected galax- 



ies but stronger than the scaling (1 -H f)^-'^^'^-'^ measured by Le 
Floc'h et al. (i2005i) to z ~ 1 for 24yum-selected galaxies. 
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